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ABSTRACT 

Elementary  waves  are  defined  as  scale  invariant  solutions  of  a 
nonlinear  hyperbolic  conservation  law  which  are  stationary  in  a  suit- 
able frame.  We  show  by  example  the  importance  of  next  to  leading 
order  scaling  data,  which  breaks  the  scale  invariance  for  these  ele- 
mentary waves.  The  example  is  drawn  from  petroleum  reservoir 
modeling  and  is  of  interest  in  its  own  right. 

The  crossing  of  an  oil-water  bank  with  a  geological  discontinuity 
is  not  a  simple  jump  discontinuity,  and  defines  a  discontinuity  of  co- 
dimension  two.  The  dynamically  stable  forms  of  these  co-dimension 
two  discontinuities  define  an  elementary  wave  of  co-dimension  two. 
In  this  note,  we  determine  the  possible  forms  for  this  elementary 
wave.  The  solution  which  appears  to  be  generic  is  not  scale  invariant. 
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1.  Elcmcntarj  Waves  .  ^r 

Elementary  waves  and  Riemann.,80tlu|ions  are  solutions  of  a  nonlinear  system 
of  conservation  laws  (      i,^y, 

u,  +  V./(«)  =  0  (1.1) 

which  are  distinguished  by  symmetry  properties.  We  take  (1.1)  to  be  an  n  x  n  sys° 
tem,  considered  in  d  space  variables.  A  Riemann  solution  is  a  solution  m  of  (1.1) 
which  is  invariant  under  the  scale  transformation  X,t  -'Si,st  ,s  >  Q  .  An  elementary 
wave  with  velocity  d  €  /?''  is  a  Riemann  solution  which  is  invariant  under  the  group 
action  jr,f-i  +  <7T,r  +  T.  Equivalentiy  these  solutions  are  transformed  into  a  sta- 
tionary solution  (m,  i>  0)  by  a  Galilean  transformation  i,  r- 1,  r  -  ;e  -  &t,  t .  The 
elementary  waves  and  Riemann  solutions  are  said  to  have  co-dimension  j  if  they  are 
independent  of  all  but  j  of  their  space  csordinatesl  For  a  more  extensive  discussion 
of  these  concepts,  see  [1,2J.  ^^'■ 

In  this  paper  we  show  by  example  the  importance  of  next  to  leading  order  scal- 
ing data,  which  break  the  scale  invariance  for  ihese' elementary  waves.  The  scale 
symmetry  breaking  solution  contains  a  length  scale.  The  example  is  important  in  its 
own  right,  because  the  geological  layers  are  an  important  factor  in  the  modeling  of 
petroleum  reserviors. 

2.  An  Example  firom  Petroleum  Reservoir  Modeling 

The  equations  for  flow  in  porous  media  are 

V  =  -  X(4)  Vp  (2.1) 

Vv  =  0  (2.2) 

s,  +  vV/(*)  =  0  ,  (2.3) 

with  the  neglect  of  capillary  diffusion  in  the  saturation  equation  (2.3).  Equation 
(2.3)  is  a  scalar  hyperbolic  equation  which  admits  discontinuities,  known  as  oil  or 
water  banks.  We  consider  a  vertical  cross  section  of  a  three  dimensional  reservoir. 
In  this  case  \  =  \{s,x,y)  has  prominent  discontinuities  across  the  boundaries  of 


distinct  geological  layers.  The  meeting  point  of  an  oil  bank  and  a  geological  layer  is 
a  discontinuity  for  (2.1)  -  (2.3)  ol  co-dimensfbn'iwo.  Since  it  should  move  with  a 
definite  velocity  (along  the  fixed  layer  boundary)  and  with  a  constant  geometrical 
form,  it  defines  an  elementary  wave  of  co-dimeniion  two. 

The  purpose  of  this  note  is  to  find  the  possible  configurations  and  velocities  of 
this  elementary  wave.  Because  of  the  assumed  scale  invariance  of  the  solution,  all 
data  is  given  with  lowest  order  scaling  behavior  only,  and  is  thus  scale  invariant 
also.  The  discontinuities  divide  the  plane  iQto  four  sectors,  bounded  by  rays,  while 
the  saturations  and  velocities  are  constant  in  each  such  sector.  The  elementary  wave 
is  then  defined  by  these  saturations. and  velocities  together  with  the  angles  between 

the  rays. 

o:  br 

The  equations  for  a  co-dimension,  two  ekto^ntary  wave  are  of  two  types. 
There  are  jump  conditions  across  each  discontinuity  and  there  is  a  dynamical  con- 
sistency equation,  stating  that  in  some  Galilean  frame,  the  solution  is  stationary. 

We  find  three  solutions  to  these,  equations .  One  solution  has  all  flow  parallel 
to  the  oil  water  bank,  and  one  has  all  flow  normal  to  the  oil  water  bank.  Very 
specific  angles  of  the  oil  water  bank  are  required  for  these  solutions.  The  third  ele- 
mentary wave  has  data  given  to  one  degree  higher  order.  It  violates  scale  invari- 
ance, having  next  to  leading  order  scaling  data.  The  oil  water  bank  is  tangent  to  the 
layer  at  the  crossing  point,  and  curves  away  on  either  side.  This  elementary  wave 
appears  to  be  the  generic  solution. 

3.  The  Jump  Equations 

Let  us  take  the  layer  located  on  the  x  axis,  with  the  region  >  <  0  having  abso- 
lute permeability  M  =^  1  and  the  region  y  >  0  having  absolute  permeability  1.  We 
make  the  simplifying  assumption  that  the  fractional  flow  function  /  in  (2.3)  is  the 
same  on  both  sides  of  the  layer  discontinuity.  This  is  correct  only  in  so  far  as  / 
depends  on  the  fluid  properties  (the  oil  and  water)  but  does  not  allow  for  a  change 
in  the  influence  of  the  rock  properties  on  /.  With  this  assumption,  we  can  suppose 


that  the  saturation  values  are  chosen  to  be  constant  in  eack'o|)tbe  four  sectors,  and 
at  a  value  such  that  the  wave  connecting  ihcm^iaJ2.3)  is  a  shock.  Also  the  frontal 
mobility  ratio  Mg  is  the  same  in  both  layers,  under  this  assumption.  We  follow  the 
terminology  detined  by  Fig.  1,  so  that  ifSdOdfP  a^:the  unit  normal  and  tangent  vec- 
tors to  the  bank  in  the  upper  layer  and  n  and  f  are  the  same  in  the  lower  layer.  Also 
we  can  suppose  that  the  velocities  are  constant  in  each  sector,  as  defined  by  Fig.  1. 
The  jump  equations  for  the  velocities  in  each  sector  are 

<  n,  V  -  v'  >=0=<  t,  V  -  hfy  >  (3.1) 

<  li,  V  -  v"  >  =0=  <  f ,  V  -  kjf'  >'  (3.2) 
The  jump  equations  for  the  layer  discontinuity  are 

V,  =  Mv,  ,  v/  =  Mv/  (3.4) 

From  (3.1)  we  have  v  -  v'  parallel  to  t  and  v  -  My  parallel  to  n  and  so 

<  V- v'.v- A#„v'>  =  0.  (3.5) 
Similarly 

<  V- v',v-A/^v'>  =  0.  (3.6) 
Let  M  denote  the  2  x  2  matrix 


(M  0) 
X  =  [o   ij 


Then  v  =  Mv  and  v'  =  Mv'  so  that 

<  V  -  v'.  mKv  -  My)  >  =  0  .  (3.7) 

Taking  linear  combinations  of  (3.S)  and  (3.7)  yields  the  following  form  of  the  jump 
relations: 

(v,  -  v/)  (v,  -  A/„v/)  =  0.  (3.8a) 

(v,  -  V,')  (v^  -  M„v/)  =  0.  (3.8b) 
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There  are  four  solutions  to  (3.8),  two  oi  which  satisfy  the  dynamical  consistency 
equations  and  thus  define  elementary  VAvet-.  "''•'  £°^- 

4.  The  Tangential  Flow  and  ihz  Ncraci  Flvw^lemenUry  Wavea 

4.1.  M^v'  =  V  ,  V  *  0 

From  (3.1),  v  -  v'  is  parallel  to  the  bank.  Since  v  and  v'  are  parallel  to  each 
other,  V  and  v'  are  parallel  to  the  bi->k.  In  this  case  the  elementary  wave  is  station- 
ary and  the  dynamical  consistency  enuation  is  satisfied. 

4.2.  v'  =  V  ,  V  #  0 

From  (3.1),  v  -  M^v'  is  normal  to  tiiel)a!ik,  and  thus  so  are  v  and  v'.  Let 

v„  =  <  V,  n  >  ,      V,  =  <  V,  n  >  . 

r. 

Then  the  velocity  of  the  crossing  point  is  v„  /  sin  6  ,  and  the  dynamical  consistency 
equation  is 

v„  /  sin  6  =  v„  /  sin  e  , 

or 

v„  sin  9  =  v„  sin  6  .  (4.1) 

Because  v  is  normal  to  the  bank,  we  have 

" "  imT  '   '""' " """ ' 

and  analogous  v  relations.   Substitution  in  (4.1)  yields 


^  =  '«H 


or 


since  V,  =  A#  v^ .  Thus  rarz-v  •'^)  AtjIo» 

■  i  •       ' 
which  yields  the  equation 


JL. 


4.3.  Vj.  =  v/  ,      V,  =  W„v/  ,      V  9t  0  . 

In    this    case    v  -  v'  =^  0,    so    v  -  v'    is    parallel    to    t    and    consequently 

t  =  e  ,  n  =  -«,.  where  *,  and  *,,  are  the  unit  basis  vectors  in  R^.  The  bank  and  the 

layer  are  parallel.  The  bank  moves  away  froBn^the  layer  dynamically,  so  this  case 

does  not  define  an  elementary  wave. 

■rtf  .taisb 

4.4.  V,  =  v/  ,      Vy  =  Af„v/  ,      V  =)fc  0  . 

In  this  case  t  is  parallel  to  v  -  v'  and  conseque^tljf* f  ^V^;,  n  =  *,.  The  bank  and 
layer  are  normal.  This  case  is  not  dynamically  stable  and  does  not  define  an  ele- 
mentary wave. 

5.  The  ElcmenUry  Wave  which  Breaks  Scale  Invariance 

We  have  found  two  solutions  for  the  elementary  wave  equations  which  have 
V  ^  0.  All  other  solutions  must  have  v  =  0  at  the  crossing  point.  We  assume  v  =  0 
at  the  crossing  point  and  expand  v  in  a  Taylor's  series  about  the  crossing  point, 

dv    ,        dv 

and  similarly  for  v'  and  the  velocities  in  the  other  sectors.  Also  we  suppose  that  the 
normal  to  the  bank  is  taken  to  first  order  and  thus  the  bank  is  represented  as  an  arc 
of  a  circle,  near  the  layer  boundary. 

With  V  and  the  normals  taken  to  first  order  near  the  cross  point,  consider  the 
jump  and  consistency  equations  to  the  same  order.  In  general  a  solution  of  these 
equations  will  lead  to  a  dynamical  evolution  of  the  crossing  angle  and  a  Hxed 


crossing  position.  Such  a  solution  is  not  (by  definition)  an  elementary  wave.  There 
is  one  exception,  however,  which  occurs  wjien  the  bank  crosses  the  layer  boundary 
tangentially  (with  first  order  contact) .   We  examine  this  case  in  more  detail. 

To  first  order,  the  velocity  along  the  layer  is 

av 


v  =  xv,  =  x^^. 


with  similar  expressions  for  v',  v,  and  v'.  The  jump  conditions  can  be  analysed  as 
before  and  have  four  solutions.  Of  these,  only  the  third  of  the  above  solutions  is 
compatible  with  the  zeroth  order  data,  namely  n  =  -e^,  t  =  e^.  Thus 

dv.  dv'  dv^,        dv'  ay 

dx  "    dx  dx  dx  dx 

The  crossing  point  velocity  is  determined  by  the  radius  of  curvature  of  the 

bank.   Thus  dynamical  consistency  is  a  condition  on  — ,  where  r  and  r  are  the  radii 

f 

of  curvature  above  and  below  the  back  respectively. 

In  the  notation  defined  by  Figure  2,  the  dynamical  consistency  equation  is 
dx  =  dx.  Here  we  write 

dy  =  v^  dt  =^  dx  Vy^  dt  . 

But  dx=rd^  up  to  higher  order  terms  and 

dy  =  r(l  -  cosJ«|»)  =  r  ^^^  . 


Equating  these  expressions  for  dy,  we  obtain 

2 
and  so 


Vj,^  ir  =  I  <i<J> 


Because  n  =  -e  ,  we  have 


dx  =  2r  Vy^  dt 


V  ^    V  ^    V  ^    V 
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and  so  dx  =  dx  implies  r  =  f.  This  is  the  dynamical  consistency  equation.  The  equa- 
tion 

dx      ^ 
defines  the  velocity  of  the  crossing  point. 
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7.  Figure  Captions  ,y 

Figure  1.  A  bank  crossing  a  layer  discontinuity.  This  figure  defines  the  velocities  v, 
v',  V,  v',  which  are  taken  to  be  constant  in  each  sector  defined  by  the  crossing  bank 
and  layer  discontinuity.  Also  the  unit  normal  and  tangential  vectors  n,  t  and  n,  i  are 
defined  in  this  figure. 

Figure  2.  A  bank  crossing  a  layer  discontinuity  with  first  order  contact.  This  figure 
serves  to  define  dx,  dx,  d^  and  dy. 
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